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TIP  VORTICES  OF  WINGS  IN  SUBSONIC  AND  TRANSONIC 
FLOWS  :  A  NUMERICAL  SIMULATION* 

G.  R.  SRINIVASAN  AND  W.  J.  McCROSKEYf 

JAI  Associates  Inc,  Mountain  View,  California  94042,  U.  S.  A. 
fU.S.  Army  Aeroflightdynamics  Directorate  -  AVSCOM 
NASA  Ames  Research  Center,  Moffett  Field,  California  94035,  U.  S.  A. 

ABSTRACT  Thin  layer  Navier-Stokes  and  Euler  equations  are  numerically 
solved  using  a  multi-block  zonal  approach  to  simulate  the  formation  and  roll-up 
of  tip  vortices  of  wings  in  subsonic  and  transonic  flows.  Several  wing  planforms 
have  been  considered  to  examine  the  influence  of  tip-cap  shape,  planform  geom¬ 
etry  and  free  stream  Mach  number  on  the  formation  process.  A  good  definition 
of  the  formation  and  qualitative  roll-up  of  tip  vortices  has  been  achieved. 

1.  Introduction 

The  phenomenon  of  vortex  formation  and  its  subsequent  roll-up  behind  wings  has  long 
attracted  attention  because  of  the  potential  hazard  to  aircraft  that  encounter  them  in  the 
flight  and  because  of  the  noise  and  vibration  problems  that  they  cause  for  helicopter  rotors. 
Although  the  formation  of  the  tip-vortex  is  a  viscous  phenomenon,  recently  Rizzi  et  al  [l] 
have  used  Euler  equations  to  numerically  simulate  the  vortex  formation.  Their  reasoning 
for  this  success  is  that  the  artificial  viscosity,  used  in  the  numerical  scheme  to  control  the 
nonlinear  stability,  played  the  role  of  the  natural  viscosity.  More  recently,  Mansour  (2j 
used  the  thin-layer  Navier-stokes  equations  in  conjunction  with  an  elliptic  grid  generation 
scheme  to  simulate  the  flow  field  and  the  tip-vortex  on  a  low-aspect-ratio  wing.  In  the 
present  work  an  alternate  but  more  efficient  method  of  multi-block  zonal  approach  j3,4) 
is  used  for  solving  the  thin-layer  Navier-Stokes/Euler  equations  for  isolated  wings  with  a 
view  to  simulate  the  tip-vortex  formation  in  subsonic  and  transonic  flows. 

2.  Governing  Equations  and  Numerical  Scheme 

At  the  Reynolds  numbers  of  interest  in  the  present  calculations,  the  flow  on  the  wing 
can  be  considered  to  be  mostly  inviscid  except  in  a  small  region  near  the  body  surface 
where  the  viscous  effects  are  important.  Accordingly,  the  flow  field  is  solved  using  a  multi¬ 
block  zonal  approach  developed  at  NASA  Ames  [3,4]  wherein  the  inviscid  outer  blocks 
are  solved  using  Euler  equations  and  the  inner  viscous  blocks  are  solved  using  thin-layer 
Navier-Stokes  equations.  The  governing  equations  are  generally  nondimensionalized  by 
free-stream  quantities  and  are  transformed  to  the  computational  domain  so  as  to  preserve 
the  strong  conservation  law  form  of  the  equations. 
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Typically  each  of  the  computational  zone  or  block  has  over  40,000  grid  points.  The 
grid  topology  consists  of  a  H-H  grid  generated  by  the  parabolic  grid  solver  of  Edwards  [5] 
and  the  outer  grid  had  the  dimensions  of  10  x  4  x  10  for  the  wing  C  and  10  x  10  x  10  for 
the  ONERA  wing.  The  above  dimensions  are  based  on  the  root  chord  of  the  wing.  Each 
viscous  zone  typically  had  61  x  27  x  25  grid  points  and  was  clustered  at  the  leading  edge, 
tip  region  and  in  the  normal  direction. 

The  numerical  scheme  used  a  diagonalized  ADI  algorithm  due  to  Pulliam  and  Chaussee 
[6]  with  several  additional  features  incorporated  to  improve  the  efficiency  of  the  code.  The 
details  of  this  as  well  as  boundary  conditions  and  data  management  are  discussed  in  detail 
in  Refs.  [3-4,7].  With  this  algorithm,  the  CPU  time  for  each  steady  state  solution  (with 
over  160, OCX)  grid  points  in  the  flow  field)  is  typically  under  3  hours  on  NASA  Ames 
CRAY  XMP  machine.  This  is  an  order  of  magnitude  savings  in  computational  time  over 
Mansour’s  code  [2]. 

3.  Results 

Although  several  wing  planforms  are  considered  in  this  study  of  tip  vortex  simulation, 
only  selected  numerical  results  are  presented  here  for  the  lack  of  space.  The  details  of 
these  test  cases  as  well  as  the  results  are  described  elsewhere  [8], 

Typical  results  in  transonic  flow  cases  are  described  here.  Figure  1  shows  represen¬ 
tative  pressure  distributions  at  two  spanwise  stations  for  Wing  C,  which  is  a  small  aspect 
ratio  wing  in  a  flow  of  Mach  number  =  0.9  and  at  a  =  5  degrees  angle  of  attack.  The 
numerical  results  are  compared  with  the  experimental  data  of  Keener  [9].  The  results  show 
good  agreement  with  the  experiments  over  parts  of  the  wing  which  do  not  have  massive 
shock-induced  separation.  Discrepancies  in  the  separated  regions  may  be  due  to  inade¬ 
quate  spanwise  grid  resolution,  the  thin  layer  assumption,  and/or  the  simple  turbulence 
model  used. 

Figure  2  shows  a  farfield  view  of  the  tip-vortex  formation  and  roll-up  for  this  wing. 
Fluid  tracer  particles  were  released  at  varous  heights  from  the  surface  and  at  several 
locations  along  the  span  to  generate  ti.is  pic  ture.  Particles  released  near  the  tip  region 
as  well  as  the  ones  in  the  separated  region  contribute  to  the  tip- vortex  formation.  The 
fluid  particles  from  the  lower  surface  (high  pressure  region)  seem  to  go  around  the  tip 
and  get  braided  into  the  particles  from  the  upper  surface  near  the  tip  region  and  from 
particles  lifting  up  from  the  separated  region.  The  trend  of  these  particle  paths  seem  to 
first  cluster  and  then  roll-up  and  move  in-board  before  leaving  the  wing  surface.  It  should 
be  pointed  out,  however,  that  the  outer  grids  are  relatively  coarse  compared  to  the  viscous 
grids  near  the  wing  surface.  Therefore  the  sharp  gradients  of  the  tip-vortex  decrease  in  the 
downstream  wake  [10]  and  the  accompanying  vorticity  contours  at  different  downstream 
locations  shows  that  the  vortex  is  getting  diffused. 

Figure  3  shows  a  farfield  view  of  the  tip  vortex  for  the  ONERA  swept  tip  blade  at 
Moo  -  0-85  and  a  =  5  degrees.  Here  again,  the  fluid  particles  from  the  lower  surface  move 
around  the  tip  to  the  upper  surface  and  get  braided  into  the  particles  here  and  lift-off  from 
the  surface  and  roll  in-board  in  the  down  stream  wake.  In  this  as  well  as  the  previous 
case,  the  tip  vortex  stays  distinctly  above  the  wake  vortex  sheet.  The  vorticity  contours 
shown  in  this  picture  at  several  X-locations  in  the  wake  indicate  the  distorting  shapes  and 


diminishing  magnitudes  of  these  due  to  the  grid. 

In  conclusion,  a  reasonably  good  definition  of  the  formation  and  roll-up  of  the  tip  vor¬ 
tex  is  simulated  numerically  using  a  multi-block  zonal  approach  for  solving  Euler /Navier- 
Stokes  equations. 
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Fig.  1.  Typical  surface  pressure  distributions  for  Wing  C.  M ^  =  0.9.  a  -  5°.  and 
Re  =  6.8  million. 


Fig.  2.  Farfield  view  of  the  tip  vortex  for  Wing  C.  The  vorticity  contours  at  three  x-locations  in 
the  wake  region  show  the  magnitudes  and  the  shapes  of  the  tip  vortex. 
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|  Fig.  3.  Tip  vortex  for  the  ONERA  swept-tip  wing.  The  vorticity  contours  show  the  magnitudes 

[  and  the  shapes  of  the  tip  vortex  at  several  x-stations. 
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